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'- Abstract 

o 

»v^ In this paper, all the known classical solutions of plane perfect plasticity system under Saint 
Venant - Tresca - von Mises yield criterion are associated with some group of point symmetries. The 

(-«) ^ equations of slip-line families for all solutions are constructed, which permits to determine explicitly 

^^ boundaries of plastic areas. 

^^ It is shown, how one can determine the compatible velocity solution for known stresses, consid- 

^^ ering symmetries. Some invariant solutions of velocities for Prandtl stresses are constructed. The 

C^ mechanical sense of obtained velocity fields is discussed. 
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Let us consider the well known system of two-dimensional perfect plasticity [7] : 



^ |"-2fc('?!cos2r+?^sin2rUo, 

'^ ox \ ox ay J 

^ ?^-2fcf?^sin2r-?^cos2rUo, 

— Oy \0x oy J 



which follows from equilibrium equations in the absence of body forces: 



O da^^ dr^ ^ ^ dr^ + ^ ^ (2) 

^_J dx dy ' dx dy 

^^ and yield Saint Venant - Tresca - von Mises criterion 

X {a,-ayf+iT^y^Ae (3) 

C^ through the change of variables (due to Levi) : 

ax = (T — ksin 29'^, ay = cr + fc sin 29'^, r^y = k cos 29'^. (4) 

Here ax, ay are the normal components and Txy the tangential component of stress in a rectangular 
coordinates (x, y), 9'^ + -k /A is a slope angle of the maximum principal stress relative to the axis Ox, a is 
the mean compressive stress (or hydrostatic pressure) , and fc is a constant of the material. 

To construct the velocity field, which is consistent with solution of (II]) of the form a — ao{x,y), 
9'^ = 9q{x, y) one needs to solve the following system of linear equations: 

du dv\ . ^^ { du dv\ 

— -f- — sm26lo + K- - 7^ cos 26*0 = 0, 

oy ox) \dx dyj 

du dv 

dx dy 
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where u{x, y), v{x, y) are the components of velocity vector. 

System (fTl) describes the stress state of material, which is being plastically deformed and has been 
studied over one hundred years. The main contribution in determining closed-form solutions was made 
by L. Prandtl and A. Nadai at the beginning of the 20th century. These solutions are now known as 
classical ones and serve as a good approximations for the real mechanical processes. 

The symmetry method of solving differential equations was applied to system (fTj) for the first time in 
[T] . In [2] the main part of symmetries for system ([2]) , ([3]) was calculated and some known solutions were 
associated with correspondent transformation groups. Finally, the complete Lie algebra of all admissible 
point symmetries of M was determined in [19j and all conservation laws as well Lie-Backlund symmetries 
were constructed. Moreover, in the series of papers '20], [22], |24j point transformation groups were used 
to deform some known solutions and in j21j conservation laws were applied to solve the main boundary 
problems. 

Let us note, that a group analysis of the differential equations, is a semi-inverse method. It means that 
firstly some solution should be found, then one can determine the corresponding boundary conditions. It 
allows to solve statically determined problems for (U]), ([5]), i.e. when boundary conditions involve only the 
stresses, sufficient to permit a determination of the plastic region and state of stress without considering 
the velocities. The velocities can be calculated afterwards. 

The interest to the system of plane plasticity has been recently renewed. In [S] the known symmetries, 
admitted by (fl]), ([5]) were used to determine some solutions in the form of a propagation wave. In |10| 
the complete Lie algebra of symmetries for (fTl) , ^ is calculated. Two new operators of symmetries are 
determined. 

In work "S] the plastic stress states of the round sphere are analyzed from the point of view of 
orthogonal equiareal patterns, and an interesting relation between slip-line fields on the sphere and sine- 
Gordon equation is shown. 

The main goal of the present work is, on the one hand, to re-consider the relation between the 
known classical solutions and the group of point transformations (symmetries) , because some of results 
were omitted previously or were not published. On the other hand, we analyze velocities, which are 
compatible with stresses, defined by (II]) from the symmetry point of view. 

The algorithm of construction of invariant solutions is well known. An interested reader can find 
a numerous bibliography on this matter (see, for example, |16j). The main steps are: 1) for the given 
system define a Lie algebra of admissible operators, generating a group of point continuos transformations 
(symmetries); 2) determine the optimal system on non-similar subalgebras; 3) construct a set of invariants 
for representatives of subalgebras and if possible, find out the form of invariant solution; 4) substituting 
the invariant solution form into the original system, the so-called factor-system is obtained, which is 
in our case the system of ordinary differential equations; 5) solving a factor-system one can obtain the 
invariant solution. 

The paper is organized as follows. In Section [T] we consider some known solutions, obtained through 
heuristic methods and their mechanical interpretation. Section [2] is devoted to correspondence between 
symmetries and known classic closed-form solutions. Finally, the construction of velocity field, related to 
stress solution is discussed in Section |3J 

The following information will be useful. System M is a hyperbolic one. Two families of characteristic 
lines are defined by two equations: 

^^tan0^ ^ = _cot0^ (6) 

ax ax 

and corresponding Riemann invariants, which are constant along characteristics, look like 

In the theory of plane perfect plasticity, characteristic lines coincide with so-called slip lines: curves 
whose directions at every point coincide with those of the maximum shear strain-rate. Equations ([6| 
define slip-line field. As indicated in 4 , the field of slip- lines is the fundamental unknown element to be 
determined. 

In polar coordinates {r,f), system (pi), (Is]) has the form 

9cr,. 1 dTr^ (Jr - <7ip _ „ 

or r aip r 

dTr^ 1 dcr^ 2 

or r oip r 

(fT, - (7^)2 + 4t2^ = 4fc^ 



where (Jr, o^ are the radial and angular components and r^^ the tangential component of stress. 
Introducing in a similar way 

CTr = cr — fc sin 26'^, G^ — a ^k sin 2(P ^ Tr^p = fc cos 20^ 

into above system, one can obtain: 



r — 2fc r ^— cos 26*^ + ^— sin 261^ ) = 2fc sin 26*^ 

or \ or dip 

f)fjP / rifjP f)f)P 

2fc r ^^ sin 26IP - ^— cos 20P ) = -2k cos 26*^, 

093 \ or dip 



(8) 



where aP — (cr,. + (Ty)/2 = a and 6*^ is an angle between radio and slope of the slip-line. 
Two families of characteristic lines are defined by following equations: 

fir or 

= -rtan9P, —^rcoteP. (9) 

dip dip 

From system (Is]), using Q and 6 = 6'^ = 6p + ip, we have: 

r^ - 2fc r V cos 2(r - ip) + ^ sm2{e^ ~ ip)] = 
dr \ dr dip J 

^ - 2fc r V sin 2{0^ ~ ip) - ^ cos 2(r - ^) ^ 0. 
dip \ or dip J 

1 Some notes about known solutions 

1.1 Revuzhenko solution 

Revuzhenko in his paper |18j considered the so called limiting equilibrium equations: 

I da dip _ I da dcj) _ 
2kdXi~dXi~ ' 2kdX2^dX2~ ' ^ ' 

dy / , Tr\ dx dy / , 7r\ 9x 

= tan ( ^~" 



5Ai 



( , ■K\ Ox Oy / , 'K\ Ox , „, 

where </) = ^"^ + 7r/4 and Ai = 77, A2 = ^ are characteristic coordinates (It]). Let us take 2fc = 1, then the 
first two equations (11) give 

a-e' = F{0,a + e' = G{i^), (13) 



and equations (12) obtain the following form 



f^=tanrf^, |.-cot^=|, (14) 

d-q drj dt, 04 

where 0'^ = [G — F)/2, G{ri) and F(0 are arbitrary functions and not identical constants. 

Solving linear equations ( |l4| ), closely related to ([6|, one can find parametric equations x — a; (^,77), 
y — y(^, 77) of characteristic lines in Cartesian coordinates. If the slip-line field is known, then one can 



calculate functions a{^,r]), 0{^,r]) at the point {x,y) using (13) 



In polar coordinates, system (14) takes the form 



^ dr ,^^ ^dif Idr ,^^ ^dp ,^^, 

--=cot(r-^)^, -- = -tan(r-^)^, (15) 

and is similar to ^. 

Eliminating r from the above system, we obtain one equation for the function u = tanOP; 

^'111 G' d IF' d ^ ,_, 

ln|u| — - — I —-U — 0. (16) 



d^dr] ' ' 2 d^u 2 drj 



Its solution can be determined by separating variables and has the following form (a, b,c — const. 



u^ = tan^ I 



^^(0 



Then, for the simplicity, one can take 



(17) 



tan-' 0P 



e 



, G = 2(772-71/8), F^2ie+7^/S). 



Regressing to variable (p = 9'^ — 9^ = {G — F)/2 — 9^ we have: 

(/3 = 77^ — ^2 — 7r/4 =F arctan -. 



(18) 



(19) 



The sign =F follows from the quadratic type of (17) 
Taking ta.n9 



±77/^ and integrating (15), which take the form: 
Idr 



C d(f 1 dr 
r drj r] dr]' r d^ 



■q_dip_ 



±2. 



one can determine function r 

r = e 

Solution for cr, 6 is given by functions G, F: 

G + F ^2 2 
-'i^ + V 



^^"V^ 



G-F 



v'-i 



2 ■» ■ ' ' 2 

From the necessary condition for an envelope of a family of curves 



dip dp 

equations (20), we have tttt;— — 0, or 

at, o-q 



4' 



(20) 
(21) 

(22) 
= 0, taking into account 



2eT 



$2 + 772 



0, 277 T 



$2+7^2 



0. 



Let us note, that two solutions can be obtained one from the other through the change 77 o —t] (or 
^ O — Oj the same is valid for conditions for envelope lines. So there is no principal difference between 
mechanical interpretations of these two solutions. 



Let us consider solution (19), (21) with the lower sign. Adding the condition 

dip 77 



m-'^^e 



rj^ 







(23) 



and expressing 77 = ?7(<^), one can obtain the corresponding envelope for the family of characteristic curves 
^ = const. Corresponding envelope for the second family of characteristics is defined by condition 



dip ^ 

-11 = 277 + ^. = 0. 

077 ^^ + 77-^ 



(24) 



Let us note, that the above equation through the change ^ O 77 is reduced to (23) and it defines cusp 



(return points) curve for the family ^ — const. Moreover, from relations (23) and (24) it follows that 
along an envelope, the product £,ri should be negative. 



1 Let us consider case ( p3| , when 77 e (—1/2,0), £, G (0,1/2). This case corresponds to envelope OB' 
(see Fig. IT]). In a similar way one can obtain envelope OB", taking the root of (24) and considering 

77 e(0, 1/2), Ce (-1/2,0). 

Right angle plastic area, bounded by curves OB' and OB" is, in some sense, an analogy of Nadai 
solution for the converging channel with straight line borders. Point O is the singular one. 



2 Taking the root of (23), namely r]{S,) = (—1 — a/1 — 16$'*)/(4$), one obtain the envelope AC (see Fig. 
[2]), which is a part of a spiral. Another boundary AB of this hornlike plastic area is one of the slip lines 
— l/2<$ = $o<0. The contact of curves is achieved at point A. 




Figure 1: Angle area of the first kind 





Figure 2: Hornlike area 



Figure 3: Angle area of the second kind 



3 Let us consider envelope OB" of the family -q — const. (Fig. pi), obtained from (24 1, taking ^{rf) = 
(-1 + yjl - 16?74)/(47?) and 77 e (0, 1/2), ^ e (-cx), -1/2). Another boundary OB' of the angle area is 
the slip line 77 = 1/2. 

For more details for boundary conditions see [18] . 



1.2 Nadai solution for two concentric circles 

In [13] Nadai proposed the form of dependence r^ip = h{r) = kcos29P, that from (Is]) gives: 

cos2eP = Cir-^ + C2, (J = -2kC2(p + f{r), 

where /(r) can be determined by quadrature from equation 

/' - 2k{ePy cos 26IP = 2kr-^ sin 26^^. 



(25) 



(26) 



This solution is interpreted as stresses in the area, bounded by two concentric circles r = a, r = 6, so 
that 



which gives 



Ci = -2.^^, C, = « 



62 -a2 



62 -a2 



(27) 



For simplicity, let us take a == 1, 6 = a/2, then, taking into account that 

1 /2-r2 
a/2 V r2 - 1 

from the first equation for slip- lines (|9| one can obtain the relation for the first family (1 ^ r ^ \/2) 

—pz arctan — — — , — arcsin(2r^ — 3) = v2o9 + Kt, 

V2 2y2V(r2-l)(2-r2) ^ ; v- i, 

which are parts of epicycloids. 

The second equation gives the equation of the second family of slip-lines (hypocycloids) : 

37^2 _ ^ 

v2 arctan — — — , — arcsin(2r^ — 3) = 2v2v' + Ko- 

Let us note, that concentric circles are envelopes for corresponding families of epi- and hypocycloids 
(see Fig. l4|. 




Figure 4: Epi- and hypocycloids 



2 Group analysis and invariant solutions for stresses 

It is known [19] that system (IT]) admits an infinite algebra of generalized (highest) symmetries. The basis 
of Lie algebra L„g of point transformations is formed with the following operators: 

d d ^^ d d d ^^ d 

ax oy ox ay at) ocr 



d d 

Xa ^ii{x,y,a,e)— +£,2{x,y,a,e) — 



'''i^-^w 



where 



X5 = e(.,^)|:+.(.,^)|, 



^1 — X cos 26 + y sin 29 + y — , £,2 = x sin 29 — y cos 26 ~ x — , 



and (f , 77) is an arbitrary solution of the linear system 



|-2fc(|cos2.-,|sm2.)^0, 
|-2fc(|sin2^-|^cos2^)=0. 



(28) 



(29) 



obtained from (fTj) applying hodograph transformations of the form x = x{<t, 8), y = y{a, 9) and assuming 
that Jacobian of transformation is not equal to zero. Operator X^ forms infinite-dimensional subalgebra 



of Lcro and we consider it later (see subsection 2.6). 

In polar coordinates above operators (except X^) take the following form; 



Xf 



r- X" 
or 



A I A XP- — 



d 



XP = rcos2{9~ip)— + (sm2ie - ip) - y) 
or \ k/ 



a d 



cr\ d^ _ AuaJ^ 

dip da k do 



Non-zero commutators of ([28| are as follows 

[X2,Xi] = -4fcX3, [X-i,Xi] - -X2/k. 



(30) 



(31) 



From the symmetry point of view, it is necessary to construct the so called optimal system of non 
similar one-dimensional subalgebras |16) in order to define different exact solutions. The optimal system 
of one-dimensional subalgebras for finite part of L^e is the following one (a € M): 



91 = (^4 + aXi) = (X4 + ardr) , 63 = {X3 + aXi) = {d^ + ardr) , 

63 = (X2 + aXi) = {d^ +dg + ardr) , 

ei±) = (X3 ±X2+ aX{) - {d, ± [d^ + de) + ardr) , Og - {X{) = (rdr) , 

a .... . 

here we use dt ^ ^r- for the simplicity. Non similar subalgebras correspond to different values of a. 

i-i . . . . 

Let us note, that system (nj) is invariant with respect to discrete symmetries: 

X — )■ — x, a — > —a, 9 — > —6; 
y-^-y, a ~> -a, 6 ^ -0, 



(32) 



and systems ([8| and ( 10 ) admit transformations 

If-)' -If, a - 



> — r 



(33) 



respectively. 

The transformation group, corresponding to operator X4 has been calculated and called " quasi-scale" 
transformation in ^22j. 

2.1 01 : quasi-scale transformation 



From the group analysis point of view, equation ( 16 1 admits the following symmetries 



F 



G 



Yi - -^di + ^c^r, + udu, Y2 = —d^, Y3 = —d,j, 
and there is no extension of the group for any specific form of G and F. For functions G, F, taken as in 



(18 1, above operators look Hke 

Yi - -^d^ + r,d^ + 2udu, Y2 = r'a^. Fa = r^'^dr,. 
Invariant solution corresponding to operator Yi has the form u = rf' f{z), where z — ^r]. Substituting 



this form into ( 16 1 one can obtain its general solution in terms of Bessel functions. In particular, / = ±l/z 



is one of the solutions, therefore solution ( 18 1 u = ±77/^ is the particular one and is invariant with respect 
to Yi. 



Taking into account solutions (19), (21), (22), we can express operator Yi in terms of r, </?, 9'^, a. 



Namely, for the higher sing solution tan^^ — 77/^ we have: 



2udu = -2^9? + 2779^, 
?7^ 



d^ = ri2rj- 
d„ = r f 2^ - 



ae+v') 
e 

V (C^ + 1]'^) 



-2e- 
277- 



V 






d^ + 2^d„ - 2^dg 



^2 ^ ,^2 



d^ + 2r^d„ + 2j^de 



Taking into account, that 



COS261P 



1 - tan^ eP ^2 



1 + tan^ OP ^2 ^ ^2 



sin2r 



2tan6lP _ 2^r/ 
1 + tan^ 0P ~ ^2 _|_ j^2 ■ 



we finally obtain 
1 



Yi = -r cos 29Pdr + [2a - sin 20P] d^ + 26" da + 2ade 



:da = -XP 



-XP. 
2 3 



For the solution with lower sign, the operator Yi is transformed by analogy. 

That is why one can associate Revuzhenko solution with sub algebra 0i, when a = 0. The case when 
a 7^ will be considered in the future work. 

2.2 ©2: (T - translation and scale 

Let us take a = — l/(2fcc) ^ 0, then invariants of 02 are: Ii — ip, I2 ~ 0, I3 = a — 1/alnr and solution 
has the form 

a = -2kc\nr + f{ip), B = B{ip). (34) 



Substituting the above form into (10), one can obtain the following system of ordinary differential equa- 
tions: 



6i'sin2(6'-v') = -c, f ^ 2kccot2{9 - ip). 



(35) 



Integrating the first equation one can find the expression for — ip, then integrating the second equation, 
we obtain (q S M) 



(T = —2kc In r — fcc In [c + sin 2(9 ~ (/?)] + const, 



= (p — arctan 



v^ 



tan 



ip + arctan tanh 



VT 



ClJ 



C2 



C2>1, 



c2<l. 



(36) 

(37) 

(38) 



Slip lines equations are the following ones: 
r = ki 



^c + sin 2(6*- ip)' 



^ko 



J/c 



^/c + sm2{9-ip) 



h e 



and tp = (p{9) is defined by (37) or (38). 

The value of the constant c is important for mechanical interpretation of stresses. 

In the case c > 1, the families of slip lines have the envelopes, and they represent the well-known 
Nadai solution |3j, [T^] and describes the flow of plastic material through the wedge-shaped converging 
channel (total angle 2a) 



n 
"+4 



Vc 



1 



arctan \l , a G f 0, — ) 



((1 2 
9q ' 



Two straight lines (p = ±a are boundaries of the channel. For more details about this solution see [24] . 
In the case c^ < 1, there is not any envelope for slip line families. But solution (36 1, (38 1 can still be 

,2) (1,2)\ 1 

(Po ' J along 

straight-line borders ip = ip^ — const., because 9p — const, along the straight line ip = const. The 

(1 2) (1 2) J l 

relation between 0q ' and lPq ia given by ( 38 ) . 

Unfortunately, a new partially invariant solution for the angle 9 announced in jlOj coincides with 
the above well-known Nadai solution for the channel. Compatible velocity solution can be found, for 
example, in [J]. 

In the case of c = 1, from (35) and after applying (33), we have a singular solution of the form 



V- 



-, 0- = 2fcln— + k-p. 



(39) 



which is a well known Nadai solution describing plastic state around a circular cavity of radius R, situated 
in an infinite medium loaded by uniformly distributed pressure p, with the tangential stress equal to zero. 
For more details about this solution see [22]. Let us note, that solution, obtained in [10] as partially 



invariant solution, corresponding to the subalgebra generated by operator K ~ X4, coincides with (39) 
As for the velocity solution of ((sl), compatible with Nadai solution see, for example, [53]. 
The non singular solution of ([35| is 



a ^ -2fclnr- fcln[l + sin 2(6*- ip)] + const., ip = + arctan ( 1 + j , A G M. 

The case a = does not produce any invariant solution, because the necessary condition of existence 
of invariant solution [16] is not satisfied. 

2.3 ©3: rotation and scale 

a) The case a = corresponds to the rotation subgroup {d^ + dg) and has invariants: /i = r, I2 = cr, 
I3 = 6 — If, then invariant solution is as follows 

o- = f{r), = ip + g{r). 



After substituting into ( 10 1 we have 



2fc 
r/' = ^— , rg' = cot2<?. (40) 

su\2g 



The second equation gives two solutions 



^ ^VC2r4-l 

Gr-^ — 1 



Taking the boundary condition at r = i? as 9{R) = ^ — i^ /2 we obtain C = 1/ R and co s2q — —R/r 



so C — ~1/R^. Then sin2g ^ ±\/r^ — R'^/r'^, r ^ R. Integrating the first equation of (40 1 and taking 
the boundary condition (t{R) = —p = const, we have 



a = ±fcln|^^ + y^-lj-p. 

For function 9 

1 / R^\ IT 1 i?2 

arccos — 

2 2 r 



= 'P+- arccos [-^1^] =^+ ^ - ^ arccos -^. (41) 



Finally, taking tan^ = (r^/R^ + l) / y^r*/R'^ — 1 and applying reflection (33) to (41 ) we obtain the well 
known vortex flow Nadai solution 

a = — fclntan - arccos —rr H — — P, 

TT 1 i?2 

0-(^-- + -arccos^. 
The boundary condition for the above solution is Tr^\r=R — kcos26P — kcos2g{R) — —k, cr|r=fl = 



—p. The homotopy of the above solution with Prandtl solution ( 50 ) was analyzed in [21] . The gener- 
alization of Nadai solution was obtained by Mikhlin [TT] for Tnp = q, \q\ < k, when r = R. For the 
corresponding velocity solution see [25] . 

b) In the case a ^ 0, the invariants of the operator {d^p +80 + ardr) are the following ones: 

h=z^ re""'^, /2 = a, /3 = -0 + ^ + -^ Inr 

2 Za 

and the form of the invariant solution is 

a = f{z), 0- ^ ^ + ^Inr ^ ^+ ^Inz ^ 9P ^ ^\nz. 
2 2a 2a 2a 



Substitution of this form into system ( 10 ) gives a compatibility condition a^ = — 1 . That means z = xzLiy 
and there is no real solution. 



2.4 64 : a-translation, rotation and scale 

a) Let a = 0, 64 = {da — d^p — do). The form of the invariant sohition is 

a = -A^ + fir), r = ^ + .g(r), A > 0, 



which coincides with (251, taking into account that g{r) = 9'P{r). The Nadai solution for two concentric 



circles is an invariant one with respect to operator of 64 . 

b) If a = and 64 = (da + d^p + de) is taken, the form of invariant solution is 

a = A^ + f{r), e- = ^ + g{r), A > 0, 



which is dual to solution (25), (26 1, (27) in the sense, that families of characteristics change between 
themselves, and there is no new solution. 

c) a 7^ 0. In such a case, the invariant, relating independent variables is A = re"""^. Taking into 
account, that Ad^ is a symmetry, let us consider the following form of the invariant solution 

a-A^ + /(A), r = ^ + 5(A), (43) 

where A = const. 7^ and g = 9^. Let us note, that ii A = ±2fc there is a simple- wave solution. 



considered later (see paragraph 2.5 ) 



Substituting (43) into (10) gives two conditions on functions / and g: 

, A 2k /A + cos 2g — a sin 2g 

A 2a cos 2g + (1 — a^) sin 2g ' 
, 1 A/(2/c) + cos25 - asin2g 

A 2a cos 2g + (1 — a^) sin 2(7 ' 

which are two separable equations with a closed form solution. 

For simplicity let us take a — 1. Then function g is defined implicitly by relation 

1 X f 2 cos 25 

lnA= / ——- — ^ an = loiq) + const. 

J A/{2k)+cos2g~sm2g ^ ^^' 

Relations ^ along characteristic lines take the form: 

d In r d In r 

— , — = ^ tan g, — - — = cot g. 

dip ' dip 

dlnr din A in r -1 

But — - — = — h 1, so tor the nrst family we have 

dip dip 

and the second family is defined by 

where Ci are the constants of characteristic lines and functions /i and /i are the following ones: 

2 cos 2g dg 



A (9) = 



A/ (2k) + cos 2g ~ sin 2g tan <? + 1 ' 
2 cos 2g dg 



A /{2k) + cos 2g — sin 2g cot g — 1 
In Fig. [5] one can observe two families of slip- lines and their corresponding envelopes: 

/ , n-^/2 1 ^ 1, \/2 + l\ 

r = exp 03 H = — — In — = 

^y 4 2(2 -V2) 4 V2-I; 

r = exp f (^ - - - - In 2 

which are two logarithmic spirals. To simplify the above integrals we take A — 2v2fc, then functions /i, 
1 1 and I2 can be expressed in terms of simple functions. 

It seems that the spiral solution was first mentioned by Hartmann |3j . Unfortunately, this work was 
not published, but one can find some results in |15| . In works [T], [2] the stresses were described. Here 
we give the corresponding slip-lines field. 

10 




Figure 5: Spiral solution 



2.5 ©5: scale transformation 

The form of invariant solution is 



a = aM, e = e{^). 



From the first equation of ( 10 1 one can obtain 

6 — (f = T:n/2, n G Z. 
The second equation gives 

(J = {-iy'2kip + const. = (-l)"2fc6l'^ + const. 



(44) 



In other words, there are simple wave solutions, which are well known [7 and which are called simple 
stress states. Let us note, that in Cartesian coordinates, these solutions correspond to the so-called 
similarity ones (or dimensionless), i.e. depend on y/x variable only. Unfortunately, it was not seen in 
llOj . Solution, obtained by the symmetry reduction for the representative sub-algebra of new symmetry 
Bi = —vdx + udy has the form (44). 

For example, in case a — 2k9 + const., from (fl]) function 9{x, y) is defined by the following relation: 



xcos9 + ysiii9 = $(0), 



(45) 



where ^{9) is an arbitrary function. Of course, one can take different forms of $ (as in [5], [TU]) to 
define some solution, but all these solutions correspond to one characteristic family, which is a family of 
straight-fines 9 — d: 

xcosCi -fysinCi = $(Ci), (46) 



with relation (14) along characteristic dy/dx = — cotCi. From (45), taking into account the second 
relations along characteristic dy/dx = tan 9, one can obtain linear equation 



dx 
d9 

where f{6) = $(6')/ sin 6*, with solution 



X = sm ( 



cote{x + f {9)81-0^9) =0, 



/ /' cos 9d9 + '''), T — const. 



The family of characteristics (45), (47) is orthogonal to the family (46). Equation 

x + f'i9)sm^9^0 
defines the envelope of straight-lines. All above relations are due to [5]. 



(47) 



(48) 
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The velocity fields for simple stress states have some trivial properties [23] ■ To construct these fields 
let us use components U and V of the velocity along characteristic directions (|6]): 

u = U cos 9 — V sin 6, v = U sin 9 + V cos 9. 



The velocity component along each straight line is constant. Thus, in the case of solution (45 ) U = const. 
along characteristic and for V from ([5]) we obtain equation 



dV . . dV 

— — sm 9 — -— cos 9 = 0, 

ox ay 

with general solution V{x, y) = V{Ji), where V is an arbitrary function of Ji = xcos( 
u = Ui9) cos9 ~ V{Ji) sine, v = U{e)ain9 + V{Ji) cos9. 
If the second family is a family of straight lines, then by the analogy, we have 
u = U{J2)cos9-V{9)sm9, v = U {J2) sin 9 + V (9) cos 9 , 

where J2 = xsin9 — ycos9, and U, V are arbitrary functions of their arguments. 




ysin9. Finally, 



Figure 6: Simple wave solution 

Let us note, that the form of the solution, corresponding to simple wave can be quite complicated. 
For example, one can verify that 



a = 2k9'^ + const, rcosi 



^) = Ce^ 



(49) 



is a solution of (10 1. Then, $(0) == Cexp(6''^) and from (48) one can obtain a logarithmic spiral 



as the envelope for the family rcos(Ci — Lp) = Ce'-^^ . The second family of characteristics is defined by 
the following relations: 

r cos (^ = Ce"' (sin 9" + cos 9"") + r sin 9", 

rsimp = Ce^ (sin 9'^ — cos 9'^) ~ t cos 9'^. 

In Fig. [6] one can observe the corresponding slip- line field. 

2.6 Infinite part of symmetries 

Classic Prandtl solution TT for system (l2|, (Is]) has the form 



ay mx 



m?y'^ 



ijcy_ 

k 



my 

X' 
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where /i>0, O^m^l and c are constants. This solution describes stress state of a thin plastic block of 
the height 2h, compressed between rough parallel plates and it satisfies the following boundary conditions 

Txy\y=±h = ±mk. 

Let us take c — 0, 2k — 1, h = 1 and to = 1 (corresponding to perfectly rough plates). Then in terms 
of functions a, 6 we have: 

2a = -x+ v/l-j/2, y == cos 261. (50) 

This solution is invariant with respect to sub algebra of the form {d^ — 2dx), where operator dx is of X^ 
form. 

Equations for slip-line families {9 is the parameter) look as follows: 

X = ^29 + \/l — y^ + const. , y = cos 29, 

and represent two orthogonal families of cycloids. 

Let us consider the solution, proposed in 9J in a form of propagation wave 

c^ = /(0>^ = 9{0^i = aix + a2y, 

where a^ are constants. It is easy to see, that this solution is invariant with respect to the operator 
X = a2dx — aidy, which is of X^ form. Lie algebra L5 with the basis {Xi,X2,X3,dx,dy), considered 
in [5], has the non-similar sub-algebra 9o = (93,-1-7X3), 7 € M. If ai 7^ 0, then by the rotation 
transformation, ai — )■ and X ^ dx (z 6ol7=o- That is why one can take ai = 0, a2 — 1, and the form of 
invariant solution is cr = /(y), 9 = g{y), that gives the trivial constant solution. 

Another invariant, considered in ^ is t = a + aio; -I- a2j/, which corresponds to operator X = 
2aidc! — dx — 01/02 dy. In the same way we obtain X ^ 2da — dx ^ 0ol7=-2, that produces the Prandtl 



solution ( 50 1 , so there is no new invariant solution. 



3 Velocity field 

Lie algebra of operators, admitted by (IT]), (l5| is known (see [2], [inj)- Substituting solution a = (Jq{x, y), 
9 = 9o{x,y), one can find the symmetries of ([5| only. Analyzing the structure of the optimal system 
of non-similar subalgebras, one can construct corresponding invariant solutions. As an example, let us 
consider the construction of velocity field for Prandtl solution. 



From ( 50 ) we have 

X = —2(7 — sin 29, y ~ cos 29. 

Taking into account that da = -~2dx, dg — -'2ydx + 2\J\ — y^dy, we have the following Lie algebra 
huv of point symmetries: 



^4 = 



Zi = udu -I- vdv, Z2 = -2ydx -\- 2\J\ - y'^ dy - vdu + udy, Z3 = dx, 
2arccos2/ + 2y (x — \/l — y'^ j 9^, + 2 ( — a; + \/l — y^j \J\ — y'^dy — {xu + yv)dy+ 



du, ^5 = ydu - xdy. 



-f yu + V ix — 2\l\ — y"^ j 

Zq = ua{x,y)du +vo{x,y)dv, 

where uq, vq is a solution of ([5]). As in the case of Lie algebra L^e, algebra L^v contains the infnite- 
dimensional operator Zq, which corresponds to the principle of linear superposition of solutions of linear 
equations. 

Let us note, that Z^ ^ Zq, because u — y, v — —x is the trivial solution for any form of ctq, ^q- 
Symmetry Z2 is an analog of rotation group and Z4 is an analog of quasi-scale transformation. 

The non-zero commutators for the finite part of Lie algebra L^v are the following ones: 

[Z2, Z4] = -4Z3, [Z3, Z4] = -Z2, (51) 

and for infinite part we have 

[Zi,Zt^] = —Zr-,, [Z^,Zz] = —dy, 

[Z2, Z5] = {-X + 2^1 - yA du + ydy, (52) 



[^4,^5] = (x^ -3y^ - 4xv/l - y2 + 2) 9„ + {-2^xccosy-2xy + 2y^/l^^\dy. (53) 
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3.1 Infinite subalgebra 



It is interesting to note, that relation (52) produces the well known Nadai solution for velocities [1]: 



2^1^^, V = -y, (54) 



and from (53) we have (Ci = const. 



X 



2 



3z/2-4xVl-y2 + Ci, 



/ ^ (55) 

V = 2y\/l — y^ — 2arccosy — 2xy + C2. 

Using scale transformation u ~ au, v = av, corresponding to operator Zi, one can multiply the above 
relation by the same factor a 7^ 0. And with the help of translations du^dy G Zq we can add any constant 
to the velocity components. 

In general, any combination of commutators of the form [.^2,4, [^2,4, • • ■ [-^2,4, Z^] . . .] produces solu- 
tions for velocities. Thus, commutator [Z2, [Z4, Z5]] gives 

u = 2xy — 2arccos2/ — 2j/-\/l — j/^, v — —x"^ — j/^ + 6, 
that coincides with Ivlev-Senashov [5] solution of the form 

u = xy + arcsiny — y\/l — y^ + consti, v = —{x + y )/2 + const2, (56) 

taking into account that arcsiny — tt/2 — arccosy. Let us note, that the above solution is invariant with 
respect to the operator (Z3 + Z5), and Nadai solution is the invariant one for the operator (Z3 + (9„). 



Let us give a mechanical interpretation. Nadai solution (52) satisfies linear boundary conditions for 
u on the edges h — ±1 of the plates 

u(x, ±1) — X, V = =pl, 

while u is a constant. 

If Ci = 3, C2 = TT, then solution (55) satisfies simple symmetric boundary conditions on the plates 
2/ = ±l: 

u{x, ±1) — x^, v{x, ±1) = t{'2x — tt). 



in other words, solution (55), in the sense of boundary conditions, is a generalization of Nadai solution, 
because now for the velocity u we have quadratic dependence on x, and for v the dependence is linear. 

One of the main properties of plastic deformation is indicated in [7] : the dissipation must be positive in 
plastic zone, because a plastic deformation is accompanied by irreversible energy consumption. This is a 
condition of compatibility of stress and velocity fields. The following condition ensures the non- negativity 
of the plastic dissipation of energy: 

du dv dv du 

dx dy _ dx c?j/ „ (57) 

If the above inequality is not satisfied, then boundary co ndit ions are wrong formulated from the mechan- 



ical point of view. Thus, for Nadai solution condition (57) is satisfied for any x: 2/-\/l — j/^ > 0. For 



solution ( 55 ) the plastic dissipation is non- negative if 

X > 2^1^^, (58) 

i.e. when Nadai solution component m ^ 0. 

To analyze the complete solution of (fTl), ([5]) one needs to know both slip-line field and velocity field. 
To construct streamlines for the known components of velocity vector field, we can solve equation 

^ = ^. (59) 

u V 

Thus, for Nadai solution the equation of streamlines is 



yy'l — y"^ + arcsiny + const 

y 
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and for solution ( 55 ) it looks as follows 



1x arccos y + yx — 2xy\/\ — y^ — t^x — y + 3y = const. 

The theory of plane plastic strain does not involve the time and calculated stresses do not depend on 
the rate of strain. The progress of the deformation can be expressed by monotonically varying quantity: a 
load, an angle or, as in case of Prandtl solution, one can take the height of the block (see [23]). Streamlines 
coincide with the path of material particles and they help to visualize the flow. One can observe vector 
fields and streamlines for Nadai solution and for solution ( 55 ) in Figs. It] and Is] respectively. 




Figure 7: Nadai velocity solution for Prandtl stresses 

Due to linearity of velocity equations there is a infinite number of their solutions. The following 
condition can be useful for selecting the appropriate one. It is known ^, that for the real velocity field, 
due to modified Mises criterion, the dissipation 



, du 

where e^. = -— , e„ 
dx ^ 



dv 
dy 



,27: 



U — ^x^x ~r ^y^y ' ^^xy^xyi 

dv du 



xy 



(60) 



dx dy 



are corresponding strains, should take a maximum value 



Dr 



> D. 



possible- 



Thus, for Nadai solution D = 1/v 1 ^ y^ and for solution (55) D = 2 



2/r 



y 



n/T 



therefore, when 



y 



X > 1/2 + 2-\/l — 2/2^ the dissipation of solution (55) is greater than dissipation of Nadai solution. 

Let us note, that Prandtl solution is a good approximation for experimental data only at sufficiently 
great distance from the centre [1]. The same is true for velocity field. 




Figure 8: Velocity solution (551 for Prandtl stresses 
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3.2 Invariant solutions 



Let us consider the finite part of Lie algebra L^y. Noting, that commutators (51 1 coincide with (31), 
when fc = 1, one can conclude that the optimal systems are similar: 

Gi = (Z4 + aZi) , 62 - (Z3 + aZi) , 63 = {Z2 + aZi) , 
ef ^ = (Z3 ± Z2 + aZi) , 65 = (Zi) . 

3.2.1 63: rotation symmetry 

Let us consider a = (only rotation operator) and introduce new unknown functions p{x,y), ip{x,y): 

u = psmip, V — pcosip, (61) 

then operator Z2 takes the form 



Z2 = -2yd^ + 2v/l - y^dy - d^. (62) 

Its invariants are Ii = z — x — y 1 — y'^ , I2 — tp ~^ — arcsiny, I^ — p and the form of invariant solution is 



P = f{z), "0 = 9{z) - - arcsiny. 
Substituting above relations to ([5]) and simplifying, we obtain system 

/g'+^/ sin 2.9 = 0, /'-^ cos 25 = 0. 
Taking g — l/2arccos/i(z) the above system takes the form 

h' = i-h\ r = i/2fh, 



and finally 



h — tanh(z + ci), f — C2 cosh ' (z + ci), Ci — const, 
u = C2 cosh ' {z + ci) sin I - arccostanh(z + ci) — - arccosy ) , 



3.2.2 Qa: z-translation and rotation 



Let us consider the case of 64 , when a = 0. Then operator of symmetry in terms of variables (61 ) is 



slightly different from ( 62 ) and looks as follows 

il~2y)d,, + 2^1-y^dy~d^. 
The form of invariant solution is 

P = f{z), ip = --Bicsmy + g{z), z = x-\/\-y'^-- arcsiny. 
To determine functions /, g by analogy with the above section we come to the system 

\ - sin 2.9") /' - jg' cos 2g = 0, /' cos 2.g - (^ + sin 2^ fg' = //2, 



with solution in quadratures 



g = arctan 



/ ^3 
2 — V 3 tanh ——z + const. 
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In/ = - / cos2g{z)dz. 



3.2.3 02: X - translation and scale 

The case a = gives a trivial solution u,v = const. Let us consider a 7^ 0, then invariant solution looks 
as follows 

u = e"^/(y), « = e"^g(y). (63) 

This form was proposed in 2J, but the solution was not determined. 
For functions /, g there is the system: 



Vl - y^ (/' + «5) + 2/ (.9' - a/) = 0, a/ + g' = 0. 
Eliminating /, we obtain the linear equation for g 

g" - 2a—L=g' -a^g = Q 

vi-r 

with particular solution (when a = —1/2) of the form 



The second particular solution can be found in a well-known way and is as follows 

52(2/) = — fl + Vl - y^ + y arcsin y\ . 



Finally, solution in the form ( 63 ) is 



y^ + y arcsin y 
y ' 



uix,y) = yeK--+^^) (1 - VT^y [ci + ^ (-1 + yr 

i;(x,y) = eK-^+^i^) (l - ,/l^y c^ + '^ (l + ^/l^ + y arcsin y~) 



In the case when ci 7^ 0, C2 = 0, from equation (59 1 we obtain a family of streamlines 

2 



X 



+ const. — \J\ — y2 _|_ In I ■\/l — y^ _ 1 j 



and boundary conditions u(x, ±1) = ±e~^/^, v(x, ±1) = e^^/^. 

But the case ci = 0, C2 < seems to be more interesting from mechanical point of view, because in 
such a case the boundary conditions are: u(a;,±l) — C2(7r/2 — \)c^^l'^, u(x, ±1) = ±02(1 + 7r/2)e^'^/^, 
i.e. the plates are coming close to one another with the same velocity. 



Moreover, dissipation function ( 60 1 has the form 

n V^l-y^ - 1 + yarcsiny i(_a;+yiz^) 

^ = ^^2— -T e^^ ^ > (64) 

2 (1 - v/r^^) v/r^^ 

and is non-negative for any y G [— 1, 1], so Z? ^ along the whole plastic block. 

The velocity field and a family of streamlines are shown in Fig. (l9|. Let us note, that streamlines 
look like in the experimental data, indicated in |14| . 

3.2.4 01 : quasi-scale symmetry 

In the case, when a = 0, operator Z4 in terms of variables p, '0 (|61[) and d = —1/2 arccos y takes the form 



Z4 = [-40 + 2 cos 26l(a; + sin 261)] d^ - {x + sin 26l)9e - cos(26' + 2ip)pdp+ 
+ [sin(26' + 2tjj) + x -I- sin 29] d^. 

The invariant for independent variables is z = (a; -I- sv[\29)^ — 40'^. But there are some difficulties to find 
invariants for dependent variables p, 9, so this case is a matter of the future investigation. 
Subalgebra 65 do not produce any form of invariant solution. 
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Figure 9: Scale velocity solution for Prandtl stresses 

4 Conclusions 

For the first time, Rcvuzhcnko solution is interpreted as invariant one with respect to quasi-scale transfor- 
mation. The direct construction of mechanically significant solution, invariant with respect to quasi-scale 
transformation is quite difficult, because the general solution of the corresponding factor-system of ordi- 
nary differential equations is expressed in terms of Bessel functions. Moreover, the group of symmetry is 
pointed out, which gives the Nadai solution for two concentric circles. 

All of the known classical solutions of plane perfect plasticity system are invariant with respect to 
some group of point symmetries. One can observe that for different values of the parameters, involved in 
one non-similar sub-algebra, there are different solutions from mechanical point of view. The equations 
of slip-line families for all solutions are constructed, which permits to explicitly determine boundaries of 
plastic areas. 

It is shown, how for known stresses one can determine the compatible velocity solution, considering 
symmetries. It seems there are no advantages in looking for invariants solution using symmetries of 
complete system (stresses and velocities) in comparison to traditional way exposed in this paper: firstly 
to solve the system for stresses, substitute its known solution into the system for velocities and determine 
the invariant velocity solutions. 

As one can see, the streamlines of velocity field even for the simplest Prandtl solution can be different. 
But the slip-line field is unique and it defines the region of the plastic state, bounded by envelopes of slip- 
line families (or sometimes by characteristic line). That is why, when one gives a mechanical interpretation 
of the obtained solution it is necessary to consider both fields simultaneously. 

For the Prandtl solution the well known Nadai velocity solution is compared to some other solutions. 
New velocity field in the form of exponential function is determined. It is shown, that Nadai velocity 
solution is not preferable in the whole plastic area. There arc other velocities fields, where the dissipation 
is greater. Due to infinite number of solutions, the problem of the construction of real velocities for 
Prandtl solutions is still open. 
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